arXiv:1509.05824vl [gr-qc] 19 Sep 2015 


Non-Perturbative Yang-Mills Condensate as Dark Energy 

Pietro Dona,^ Antonino Marciano /’B Yang Zhang,^ and Claudia Antolini^ 

^Department of Physics & Center for Field Theory and Particle Physics, Fudan University, 200433 Shanghai, China 
^Key Laboratory for Researches in Galaxies and Cosmology, Department of Astronomy, 

University of Science and Technology of China, Hefei, Anhui, 230026, China 
(Dated: September 22, 2015) 

Models based on Yang-Mills condensate (YMC) have been advocated in the literature and claimed 
to be successful candidates to explain dark energy. Several instantiations of this simple idea have 
been considered, the most promising of which are reviewed here. Nevertheless, results previously 
attained heavily relied on the perturbative approach to the analysis of the effective Yang-Mills 
action, which is only adequate in the asymptotically-free limit, and were extended into a regime, 
the infrared limit, in which confinement is expected. We show that if a minimum of the effective 
Lagrangian in 0 = —pafiv j 2 gxists, a YMC forms that drives the Universe toward an accelerated 
de Sitter phase. The details of the models depend weakly on the specific form of the effective Yang- 
Mills Lagrangian. Using non-perturbative techniques mutated from the functional renormalization 
group procedure, we finally show that the minimum in 9 of the effective Lagrangian exists, thus 
YMC can actually take place. The non-perturbative model has properties similar to the ones of the 
perturbative model. In the early stage of the universe, the YMC equation of state has an evolution 
that resembles the radiation component, i.e. Wy —>■ 1/3. However, in the late stage, Wy naturally 
runs to the critical state with Wy = —1, and the universe transits from matter-dominated into 
dark energy dominated stage only recently, at a redshift the value of which depends on the initial 
conditions that are chosen while solving the dynamical system. 


I. INTRODUCTION 


Observational data collected over the last two decades from Supernovae Type I a (SN la) confirmed that the Universe 
is undergoing an accelerated phase of expansion. The first evidences for such a behavior have been discovered by two 
independent collaborations, and reported in Ref. [I|, by the High-redshift Supernova Search Team, and in Ref. 0jby 
Supernova Cosmology Project Team. Analyses from combined SN la dataset (see e.g. the one reported in Ref. Q), 
cosmic microwave background radiation (CMB), for which we refer for instance to Refs. and from large scale 

structure [1,0, have further provided consolidated evidence for the current acceleration of the Universe. The source 
of late-time cosmic acceleration has been dubbed as “Dark Energy” (DE), an exhaustive theoretical characterization 
of which is still lacking. There have been many attempts so far to individuate the origin of DE, but a consensus in 
the literature has not been reached yet. 

A review of the models so far advocated within the literature of DE is beyond the purpose of this study, and we 
prefer to address reader to a sizable and rich literature that exists on this subject (see e.g. Refs. CMl)- In what 
follows, we will focus on a rather simple idea, which address the problem of DE from the perspective of condensation 
of Yang-Mills helds. Cosmic acceleration as a source of cosmological inflation has been first proposed by Yang Zhang 
in Ref. [l3l| . and then further developed b y t he same author and collaborators in the framework of current cosmic 
acceleration, as a source of DE, in Refs. 13 and [ll]j respectively in the perturbative two-loops and three-loops 
analysis of the effective action of Yang-Mills theory. To be more specific, in Ref. [3 the author has considered a 
Yang-Mills gauge boson condensate as described by the renormalization-group-improvement (RGI) action within a 
homogenous and isotropic FLRW background. Following Refs. [l6l - [T§| - the action for SU(N) Yang-Mills fields has 
been reshuffled in terms of an effective running coupling constant g = g{T), namely 
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in which stands for the background metric and n is the square of the renormalization mass-scale. From now on, 
for simplicity of notation we dehne the contraction of the field-strength tensors 
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which plays the role of an order-parameter for the YMC, and allows us to write the effective Lagrangian in a more 
compact way r = In \6/k^\. A further simplification within the studies so far developed in the literature concerns the 
use of a SU(N) gauge symmetry group, in which the number of colors N is not fixed a priori. This choice connects 
the physics under investigation to the constituents gauge-groups of the standard model of particles. 


A. The perturbative expansion: the one-loop effective action 


Moving from the two-loops perturbative (ultraviolet) expansion at large t/k^ of the effective coupling constant, which 
from the analyses reported in Refs. [1^ [l^ reads 

g-2(r) = 6oT-b2^1nT, (3) 

bo 

the analysis within Ref. [l^ has first focused on a one-loop expansion. This latter entails considering the effective 
action 


Aflf = 2^0^* In 



(4) 


The constants bo and bi are of order one, and depend on the number of colors N of the SU(N) gauge-group. Relying 
on the form of the effective Lagrangian that has been derived for the asymptotically free regime, it has been shown 
that when the minimum is attained, namely when 0 = ~ 6 oi the energy density becomes 

The equation of state then reads like the one of dark energy, namely 

p = -p. ( 6 ) 


B. Two-loops and three-loops effective action, and the issne of stability 


The issue of proving that the dark energy behavior of YMC is stable with respect to higher order loop-corrections 
has been addressed in 0 and 0. In Ref. 0 , the analysis has resorted to the two-loops effective Lagrangian, 



In 


-b 77 In 

In 

e 

+ <5 




. 


(7) 


derived within the asymptotically-free regime and then extended to the infrared confining regime, in the range of 
values for 9 in which £eff has a minimum. Again, the coefficients b and S depend on the number of color N of the 
gauge-symmetry group, and specifically read b = and 77 = 2^, with bi = The energy density and the 

pressure from this model are provided by the relations 
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At high energies, when r ^ 1, the equation of state of the YMC evolves towards the equation of state of radiation 

w = ( 10 ) 

P 3 

Approaching the infrared regime, at energies such that t < —1.8, the weak energy condition (see e.g. 0 ) is violated, 
and the equation of state parameter crosses the value w = —1. Notice however that, moving from the expression 
of £eff determined in the asymptotically-free regime, while reaching the range of values for r that ensure a dark 
energy behavior for the Universe, a divergence is encountered in the effective coupling constant 5 ^(t). Although 
divergences do not occur in the observable quantities, which are the pressure p, the energy density p, and their ratio 
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w, the occurrence of a divergence in g^{T) makes less trivial the act of relying on the extension of this procedure to 
higher-order loop expansions. Nevertheless, for the purposes of an analysis of the dark energy behavior, it has been 
shown that all the physical quantities behave like smooth functions in the range of the infrared values r G (—2.04, 53). 
Furthermore, the stability of the system accounting for interactions with matter and electromagnetic radiation, into 
which the YMC may decay, has been successfully checked. To achieve this goal, in 0 the authors have considered: 
the dynamical system provided by the first Friedmann equation in presence of matter with energy density pm and 
(electro-magnetic) radiation with energy density pr 

= —^{pPm + Pr) (11) 

the equations of motion for the energy densities component of YMC, matter and radiation that arise from conservation 
of the total energy-momentum tensor and from the decay of the YMC into matter. In comoving coordinates, in which 
H = these latter equations read 

p + 3H{p + p) = -Tp, 

Pm F Pjyi — Fp , 

Pr + m[pr + Pr) = 0 . ( 12 ) 


The decay rate F is a free parameter of the model that enters the definition of the dimensionless dynamical system 
to be solved, equation (HU and equations (HU- The specific value of the rescaled parameter 7 = ■§’ then affects the 
attractor coefficients in the analysis of the stability of the model. Consistency with Big Bang Nucleosynthesis (BBN) 
[ 23 , with the fractional density observed for dark energy (fl ~ 0.73) and with the observational constraint on the 
equivalence between the energy densities of radiation and matter at the redshift of recombination, then set the initial 
conditions of the dynamical variables. Finally, the model has been shown to posses a dark energy tracking solution, 
and the dynamical system to have a fixed point, which is stable. For 70 = 0.5, the parameter of the equation of 
state takes the asymptotic value w = —1.14, provided that in £eff one sets up (5 = 3. A different value of this latter 
parameter, for instance 6 = 7, would rather entail w = —1.18. 

A further analysis, developed along the same lines as the one reported above, has been deepened by the same authors 
in [ 1 ^, and the investigation extended to the case of the three-loops effective action. Within this latter work, the 
effective three-loops coupling constant reads 
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and the main predictions have been substantially unaffected by the improvement in the perturbative analysis. For 
different choices of 70 , which is the decay rate parameter rescaled by H evaluated at present time, the parameter of 
the equation of state at present time takes the values wq = 1.05 (given that 70 = 0.31) or wq = —1.15 (given that 
70 = 0.67). 

We emphasize that in both cases summarized above, the two-loops expansion within 0 and the three-loops 
expansion within flU | , the dark energy behavior arises from the perturbative computation of the effective Lagrangian 
in the asymptomatically-free limit. The validity of this procedure is then extended to the infra-red regime of the 
Yang-Mills theory, in order to recover a minimum in 0 for the effective Lagrangian, and derive the equation of state 
for the Universe, which entails accelerated expansion. But the occurrence of divergences in the effective action may 
shed some doubts on the validity of the conclusions for the dark energy behavior of the theory. 

Finally, we would like to point out that the gauge interactions taken into account here might not be necessarily 
considered to be the ones constituting the standard model of particle physics. It is interesting to notice that an extra 
“dark sector” might be advocated to explain the gauge group here involved. Furthermore, it might be tempting to 
identify the gauge group copies with suitable candidates for dark matter, postulating the existence of “dark copies” of 
fermions colored under the extra “dark gauge group”. The dark matter sector that is then introduced may eventually 
be connected to Mirror Standard Model theories [ 2 ^, [ 2 ^ discussed in the literature. 


II. TOWARD A NON-PERTURBATIVE INFRARED ANALYSIS 

Within previous studies 0-0 , the existence and the stability of a dark energy mechanism based on YMC has been 
investigated. The analyses have been developed moving from the ultraviolet results recovered in the literature, up 
to three-loops in the effective action for Yang-Mills SU(N) fields, and results have been then extended in 0 - 0 | to 
the infrared regime, in order to derive a physical characterization of an accelerated expansion of the Universe. An 







4 


important technical issue is that one that concerns the stability of the result at higher loops expansion, since the 
appearance of additional term in the effective action might indeed spoil stability, which totally relies on ultraviolet 
perturbative expansion. 

Furthermore, we know that the infrared regime of Yang-Mills SU(N) theories has a completely different behavior 
than the ultraviolet regime, the latter encoding asymptotic freedom, while the former showing a confinement behavior 
that depends on the number of gauge colors involved, and even on the number of colored fermions involved in the 
theory. 

Therefore, the main motivation of this study has been to prove that under mild and general assumptions, which 
are basically the existence of a minimum in 9 in the non-perturbative effective Lagrangian, a dark energy behavior 
is recovered. Then by making use of the non-perturbative techniques mutated from the functional renormalization 
group procedure, which is more adequate to be used in the confining infrared limit of the theory, it is possible to 
prove that a such a minimum indeed exists. 

At this purpose, we provided the explicit example of the SU(2) Yang-Mills action, deriving the effective Lagrangian 
for such a particular model and deepening the cases in which interactions with different forms of matter is considered. 
The procedure, which might not be completely reliable for the precise determination of the coefficients of the effective 
non-perturbative Lagrangian, is nevertheless enough to ensure that a minimum of 9 exists for £eff, and thus that a 
YMC works as reliable candidate to explain the origin of Dark Energy. In the following sections, we unfold detailed 
arguments in support of this thesis. 


A. Yang-Mills effective action from a non-perturbative approach 


Within the perturbative YMC model for dark energy that we have reviewed in the previous section, the effective 
Yang-Mills Lagrangian is that one calculated at one-loop in the seminal work by Savvidy 


23|, namely 


Cefi log 


(14) 


where b = (llA)/247r^ for SU(N). Higher-loops correction have been then deployed in order to check whether the 
substance of the physical results remains unchanged, and how the details may vary. 


Instead of keep considering higher-loops improvements of (1141) . the main purpose of this work is to start from a general 
non-perturbative form of the action, i.e. 


CeB =W{9) . (15) 

We may then proceed to determine the form of W(0) by stating some general requirements that must be fulfilled 
in order to obtain a YMC that can work to explain the origin of dark energy. The function W will be in general 
equipped also with an energy scale k in analogy with (HD for dimensional reasons. We notice indeed that W must be 
a (not necessarily analytic) function of 9 satisfying the following properties: 

1. VV(0) has a non trivial minimum in 9; 

2. yV{9) possesses a perturbative limit, which reproduces the one-loop result derived by Savvidy [^ : 

3. W{9) shows the asymptotic behavior of being at least linear in 9, which in turn is linear to the bare Yang-Mills 
action. This latter requirement can be formalized as it follows: in the ultraviolet regime 0 ^ it must hold 
the limit 


dlnW 
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In what follows we will consider for simplicity a pure SU(2) Yang-Mills theory, the gauge field of which is not coupled 
to any other fundamental matter fields. We want to stress that SU(2) Yang-Mills fields introduced here, as well 
as SU(N) Yang-Mills fields dealt with in [l3l - [l^ . despite being suitable to build a model for cosmic dark energy 
should not be necessarily identified as Standard Model gauge fields. The introduction of additional copy of SU(N) 
Yang-Mills fields might anyway allow to make contact with some Mirror Standard Model theories that have 

been advocated in the literature of dark matter. Again, such a link is not necessary for our purposes and will not be 
analyzed within this investigation of dark energy YMC. Nevertheless it might represent some interesting directions to 
be followed in forthcoming studies. 






5 


B. YMC as cosmological dark energy 

In what follows, we discuss the cosmological consequences of the requirements we demanded above for the non- 
perturbative effective action, and shed light on the behavior of the YMC for the fate of cosmological dark energy. 
We will assume a flat Friedmann-Lemaitre-Robertson-Walker (FLRW) universe, the line element of which can be cast 
respectively in terms of comoving or conformal coordinates as it follows 

ds^ = dt^ — {t)5ijdx'‘ dx^ = — Si^dx^dx^], (17) 

where the cosmological time t and the conformal time r are related through dt = adr. We will consider the simplest 
case of a universe filled only with the YMC, and assume it to be minimally coupled to the gravity. Then the effective 
action reads ^ 
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From now on, we will only denote with g the determinant of the metric TZ is the scalar Ricci curvature, and £eff 
is the effective Lagrangian of the YMC, described by Eq. m- Varying S with respect to the metric one obtains 
the Einstein equation = SttGT^i/, where the energy-momentum tensor of the YMC is given by 
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We may set up a gauge that preserves the isotropy and the homogeneity of the FLRW background. We write gauge 
fields as functions of the cosmological time t, and choose their components so to satisfy Aq = 0 and Af = SfA{t). This 
choice indeed ensures that the total energy-momentum tensor is homogeneous and isotropic. We then introduce 
the usual definition of the Yang-Mills tensor, cast in terms of the group structure constants /“^'^, namely 

- d.A^ + r^^A^Al . ( 20 ) 


For the SU(2) gauge-group to which we are specializing our analysis, the structure constants reduce to 
Furthermore, thanks to the gauge fixing we have selected above, and looking at the case of constant electric field (we 
will assume for simplicity in the following a vanishing magnetic field), we can rewrite the Yang-Mills tensor in the 
simplified form 
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( 21 ) 


This allows to express 0 in a very simple form, i.e. 6 = X]i=i each component of the energy-momentum 

tensor then rewrites as 


wro = -ln’(e) 6 l + ^w'ie)E^ 5 °^, (22) 

(“)rj = ^yv(0)S}-^yv'(0)E^s}S2. (23) 

o o 


These tensors are not yet isotropic, their values depending on the direction of the color a. Nevertheless, the total 
energy-momentum tensor isotropic, and the corresponding energy density and pressure are 

given by 


Pymc = — W(0) -I- 2>V'(0) 0, Pymc = W(0) — - W'(0) 0 . 


(24) 


^ Following the definition in ISOl . we adopt the sign convention for the metric, the Riemann tensor and Einstein equation 

respectively. 
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Consequently the equation of state (EOS) of the YMC is immediately recovered to be 


wymc 


Pymc 

Pymc 


w-lwe i-|w^ 
W - 2W' 6 » “ ~ \-2^ e' 


(25) 


It is worth discussing the mathematical properties of the EOS wymc- On one hand, if we require the Yang-Mills 
theory to condensate, then the function W must has a non trivial minimum, as we demanded in Sec. Ill Al But this 
is equivalent to require that W' vanish at some point 9q. In this point the YMC has an EOS of the cosmological 
constant with wymc = ~1- Around this critical point the YMC dark energy models can account either for an EOS 
characterized by wymc > — 1 or for an EOS characterized by u>ymc < ~li thus encoding phantom matter behavior. 
On the other hand, in the high-energy-scale regime, in which with 9 ^ one would like to recover that the YMC 
exhibits an EOS of radiation, characterized by wymc = 1/3, in analogy with the perturbative analysis [IM3. 
Within the framework of the non-perturbative action (1151) . this corresponds to the third requirement listed in Sec. 
Ill Al The effective action should then scale for 9 ^ at least like the bare Yang-Mills action, i.e. c? In W(6()/c? In 0 ~ 1. 


In the following sections, we will show in detail that the YMC evolves from the EOS with u>ymc = 1/3 to u>ymc = ~1 
while the Universe is expanding. 


C. A non-interacting YMC model 

The cosmological model we are about to analyze in this section entails three different sources for the energy- 
momentum tensor: i) the dark energy, the role of which we assume to be played by the YMC; ii) the matter, 
including both baryons and dark matter, which is dealt with as non-relativistic dust with negligible pressure; iii) the 
radiation, the component of which consists of photons and possibly other massless particles, such as neutrinos. We 
will describe the three components in terms of their EOS, without accounting for any microscopic treatment in terms 
of fundamental matter Helds. 

Since from Eq. 113 we assumed ab initio the universe to be flat, fractional densities will sum up to one, namely 
f^YMC + -I- Uj. = 1. Indeed, fractional energy densities are defined as Oymc = Pymg!P tot, = Pm/Ptot, 

fir = Pr!Ptot I and the total energy density is given by ptot = Pymc + Pm + Pr- The overall expansion of the universe 
is determined by the Friedmann equations: 

(PYMC + Pm + Pr), 

- = -:^(PYMC + 3pYMC + Pm+ Pr + ‘iPr), 

a 6 

where the dot denotes the d/dt. 

As a first preliminary investigation, we will assume that there is no interaction between the three energy components. 
The dynamical evolution of these latter is determined by their equations of motion, which in turn follow from imposing 
the conservation of each component of the energy-momentum tensor, 

Pymc + 3-(pymc + Pymc) = 0, 
a 

Pm ~t~ ^~Pm ~ d; 
a 

Pr ~t~ 3 (pr ~\~ Pr) — 0. 
a 

From Eqs. (12911 and (I30p . we can immediately obtain the standard evolutions of the matter and radiation components. 
Pm oc a~^ and pm oc a~^. A less obvious but still rather simple task is solving the evolution of the YMC component. 
Inserting (l24l) into (|28)) . we obtain the following relation, 

9 {W + 2 W ' 9 )+ A-W 9 = Q , (31) 

a 

which is in a quite compact form. This equation is integrable for any regular enough W. In particular, we can easily 
derive the result 


(28) 

(29) 

(30) 
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(32) 
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where a is a coefficient of proportionality that depends on the initial conditions. 

At very high redshift, Eq. (15^ entails an increase of the order-parameter 0 that involves the limit 6 ^ k^. Thus, 
at very-high redshift the system transits towards the ultraviolet regime. Within this limit, Eq. (1251) encodes the EOS 
parameter wymc 1/3. The YMC then starts behaving as the radiation component, as one would have expected 
since the theory is evolving towards asymptotic freedom at high energy. At small redshift, the expansion of the 
universe requires the LHS of Eq. (|32|) to vanish, which occurs for the extremal value ot 6 = 9o, and the EOS’ pa¬ 
rameter then converges towards wymc = ~1- This ensure that the YMC behaves as an effective cosmological constant. 

Finally, notice that we may proceed as in MM. and take advantage of the observational constraint on the ratio 
between the dark-energy-density and the critical energy density, in order to fix the energy-scale that appears in the 
effective Lagrangian W(0). The effective Lagrangian is no more dependent on the parameter n, and we can rescale 6 
as in the previous literature. 


D. Interacting YMC models 

In this section, we generalize the YMC dark energy model and take into account some effective interaction with dust 
matter. Nevertheless, for the sake of simplicity, at this first stage of the analysis we will disregard the interaction 
between radiation and the YMC. We will then describe the YMC dark energy and background matter interaction 
through one additional parameter Q. The equations of the conservation of energy in (E51) and (1^ should be modified 
into 


a a 

Pymc + 3-(pymc +Pymc) = — 
a a 

(33) 

Pm H” 3 Pjyi = H (y, 
a a 

(34) 

P'P 3 \Pr Pv) — O' 

(35) 
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The interaction parameter Q, in natural units, has the dimension of [energy]^ and has been introduced for 
phenomenological reasons. Its possible form will be addressed later. If the YMC transfers energy to the matter, 
for instance if the YMC decays into pairs of matter particles, we should require the parameter Q to be positive. 
Vice-versa we should require the parameter Q to be negative. 


We can then proceed to the study of the evolution of the system (l33)) - (l34l) . It is convenient to introduce the so-called 
e-folding time N = Ina, the derivative with respect to which we will denote with a prime. We will denote as x the 
dimensionless matter density x = pm/■ The system (I55|) - (IM1) then reads 

9'{W + 29W") + ^9W = -Q, (36) 

x' + Sx = +Q. (37) 


Using the above definitions, we can immediately recast the fractional energy densities of the YMC and the dust matter 


f^YMC 


-W + 2W'9 
-W + 2W'9 + x’ 


and flm 


X 

-W + 2yV'9 + x' 


(38) 


It’s useful to discuss the general properties of this dynamical system before specifying the form of the interaction 
term Q. We can seek the fixed points of the system by imposing 9' = x' = 0 in Eas. (l551) - (I571) . and then look for the 
solutions 9c and Xc of of the simplified system derived from Eqs. dSU-dSZl)- This latter reads 

A9cW{9c) = -3x„ (39) 

ixc = +Q{xc,9c). (40) 

The stability of the solutions of these differential equations, and the possible existence of attractor solutions, will be 
analyzed in the forthcoming subsections, in which we will specialize the form of the coupling Q by assuming different 
linear combinations of the energy-density components considered so far. 




1. Q OC pYMC 


In this section, we parametrize the interaction as proportional to the YMC energy-density, namely Q = apYMC = 
a(- W -I- 2yV'9). The trivial case a = 0 reduces to the free YMC dark energy model studied above. We will just 
consider the simplest case with a being a non-zero dimensionless constant. The evolution equations (I36p - (I37|) then 
recast 


0' (W' + 20W")+40W' = -ai-W + 2W'0), (41) 

x' + 3x = +a{-W + 2W'0). (42) 

When the fractional density of the YMC is sub-dominant in the universe, we expect the effect on the dust to be 
small. Only in the latest stage of expansion of the universe, where the YMC dark energy dominates its evolution, the 
effect of interaction on the dust component can become important. 

The critical point (IMll - diUl) now rewrites 


2 (2 +a) 0cW{0c) 

SXc 


aW{0c), 

-40,W'i0c). 


(43) 

(44) 


It’s easy to see that the fractional energy densities of the YMC and the EOS at this critical point do not depend on 
the details of the potential W, but might rather have a dependence on the parameter a. It is straightforward to verify 
this from the very definition of EOS calculated at the fixed point solution of Eq. (l43l) : 


f^YMC 


1 ^ yVjOc) - 20,W'{0c) ^ 3 

WYMC W{0c) - |0cW'(0c) 3 + a ■ 


(45) 


The constraint 0 < Oymc < 1 implies that a > 0. Internal consistency of the model requires the solution of the 
system to be stable under perturbations. The next step is then to require the critical point to be an attractor solution. 
In order to achieve that, we first need to compute the eigenvalues of the linearized system (l4T]) - (l42l) at the critical 
point: 

Ai = -3, (46) 

Az =- - -2 ((4 + a)Wi0cf + (2 + a)402yv;"(0j2 ^ 4(3 ^ a)0cW'i0c)W"(0c)) ■ (47) 

iW'(0c) + 20cW"(0c)) J c y j y jj 

The solution is an attractor if 

(4 -k a)W'i0c)^ + (2 -f a)40lW”(0c)^ + 4(3 -f a)0cW(0c)W'{0c) > 0 . (48) 


In general we will need a specific form of W to further discuss the nature of the critical point. 


C OC pm 

The next case to be considered hinges on an interaction of the form Q = apm = ax. The evolution equations dMll-dSTl) 
now rewrites 


0' (W + 20W)+WW = -ax, (49) 

x' -I- 3x = +ax . (50) 

If we assume a to be a non vanishing constant, we easily recover that pm oc a°‘~^. This result might then lead to 
observational inconsistencies: the evolution of the dust component conflicts with the evolution of dust in the standard 
big-bang model. We should then avoid considering such a form of interaction term in the early stage of evolution of 
the universe, at very high redshift. 

Nevertheless, even if we insisted in phenomenologically describing at small redshift the interaction between YMC 
dark matter and matter energy-density with a term of the form Q = ax, we would find only a trivial critical point 
0c = 00 and Xc = 0. Thus, we must conclude that it is impossible to obtain an attractor solution for this kind of 
system. 
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3. Q OC pYMC + p m 


As a final case, we can discuss a model where Q = a (pymc + Pm)- We limit again ourselves to the consideration of a 
being a dimensionless constant. In the later stage, when the dark energy dominates the evolution, this model reduces 
to the case we studied in Sec. HID 11 while in the dust dominated stage reduces to the case we studied in Sec. HID 21 


In analogy to the discussion developed in the previous section, if we insisted in applying this model to the description 
of the early universe, the evolution of dust would turn to be incompatible with the prediction of the standard hot 
big-bang models. Nevertheless, if we want to develop a phenomenological model for small redshift, we may elaborate 
on this case. 

The dynamical equations in (l36|) - (l^ become 


o' (W' -k 2ew") + 40W' = -a (-W -b 2>V'6i -b x) , 
x'+Zx = -ba (-W-b 2>V'6i-b x) , 

and the system admits a critical point in 

3aW = 2(a + 6)dcW', =-^0cW'(0c) ■ 

The fractional energy density and the EOS of the YMC at this critical point now read 

Dymc = — ^ 


wymc 


(51) 

(52) 

(53) 


(54) 


Notice that now the observational constraint of 0 < Dymc < 1 sets a different range of allowed values for a, namely 
0 < a < 3. 


III. A NON-PERTURBATIVE EXAMPLE: SU(2)-YMC 

We focus now to the analysis of the YMC model we have described in the previous sections, with a specific choice 
for the IR effective Lagrangian. We keep analyzing a YM theory that enjoys a SU(2) gauge group, and show at 
non-perturbative level that a YMC forms that drives accelerated expansion of the universe at small redshift. The 
Functional Renormalization Group (FRG) approach to non-Abelian gauge theories will be particularly fruitful to our 
purposes, as it allows to introduce non-perturbative methods that can be treated as much as it is possible analytically. 


A. Functional Renormalization Group 


The FRG approach is a tool developed to study interacting QFT and statistical systems in the non-perturbative 
regime, where no small coupling exists and perturbative techniques are not applicable. The method is based upon 
a Wilsonian momentum-shell wise integration of the path-integral; a mass-like regulator function Rk(j>) suppresses 
quantum fluctuations with momenta lower than an IR momentum cutoff scale k. This allows us to define a scale- 
dependent effective action, the flowing action F^, which only contains the effect of quantum fluctuations with momenta 
greater than k. By changing k we can interpolate smoothly between the microscopic action Ffc-^oo and the full quantum 
effective action Ffc^o- The scale-dependence of the flowing action is then given by the Functional Renormalization 
Group Equation (FRGE) [2^ and 


dtTk = isTr (l 


,( 2 ) 


+ R-k^ 


dfRk 


(55) 


( 2 ) 

Herein, F), denotes the second functional derivative of the flowing action with respect to the fields, and constitutes 
a matrix in field space. The Super-Trace STr includes a summation over all discrete indices and the fields, including 
a negative sign for Grassmann valued fields, i.e. fermions and Faddeev-Popov ghosts. The Super-Trace also includes 
a summation over the eigenvalues of the Laplacian in the kinetic term. The main technical advantage of the FRGE 
lies in its one-loop form, which nevertheless takes into account higher-loop effects, as it depends on the full, field- 

dependent nonperturbative regularized propagator ^F^^^ -b (see [l^)- The FRGE has been extensively applied 

to SU(N) Yang-Mills theories, for further references see |37H4l] | and for the application of the FRG to the study of 
YMG [13^. 
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B. Finding the effective Lagrangian 

Solving the equation (1551) exactly is a titanic task. Since we are mainly interested in qualitative and as much as 
possible analytic results, we will deploy some approximations. (For the reader interested in the state of art of YMC in 
the FRG framework, we refer to the work by Eichhorn, Gies and Pawlowski [d^, in which a numerical extrapolation 
between low and high momenta of full propagators were used to compute the gluon condensate.) 

First of all we will replace F^ in the RHS of equation (1551) with the bare action S^, doing so we are allowed to 
integrate both sides of the equation 

Ffc = - y £eff = - y m{ 0 ) = y *isTr (5(2) + Rk) ' dtRk = ^STrLog 

We select the bare action to be S' = | / that corresponds to the UV limit of our effective theory. We will 

fix the integration constant requiring the effective action to vanish for vanishing field strength. 

To correctly invert the regularized propagator we need to include in the action a (harmonic) gauge-fixing and the 
associated Faddeev-Popov ghosts. 


Sgf = 


2 a , 




Sgh = 


da; D^CvD^c'' 


In the Landau gauge a —>• 0, the trace over the gauge field space is restricted to the transverse sector 

isTrLog (S(2) + Rfc) = iTrt.ansv. Log {Dl^r + R^ (5^")) - ^Tr^host Log + R^ , (56) 

with operators = OScbS^'^ + fabc and in which g is the YM coupling and the barred quantities 

are made of background fields; for more details on the actions and its variations we refer to the Appendix We will 
employ the simplest possible regulator functions (mass-like cutoff) 

Rk{V) = k‘^, (57) 

in both the transversal gauge and ghost sector V , ^gh- 

We can employ an integral representation^ of the logarithm, in order to find an exactly summed expression 

r°° rlv 

Log(A) = -/ -e-^Y (58) 

A wise choice of the background allow us to perform the traces as sums over the eigenvalues of the chosen kinetic 
operators. In general the eigenvalues of the operator Dt are not known, and the only known stable covariantly 
constant background is the self-dual one already employed in the FRG context in (the key properties needed for 
this work are summarized in the Appendix |B| . The effective Lagrangian is finally recovered to be 


Zieff 


27r^ 


nOO 1 0<D 

J _ 'y ^ ^^-2gB{n+m)+k^ ^-2gB{n+m+2)+k^ _ ^-2gB{n+m+l)+k^J 


m,n—0 




( - i - 

V4sinn {gBs) 



(59) 


The “magnetic field” B is the only variable of the self-dual background, which is related to the order parameter 
through 9 = B^. The next step is to remove the constant part from the integral that gets contribution from the lowest 


^ This kind of approximation is usually called “one loop” in the FRG literature because of the similarities between the FRGE and the 
standard one-loop effective action. 

^ One should actually be more careful with the definition of the integral. A more precise formula is the following: 

Nevertheless we will use the naive representation and implicitly regularize the final expression. 
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order expansion of the sinh. We then perform a change of variable, and reshuffle (l59l) as 




^ r ^ +1-—) 

Jq s V4sinh^(s) ds^y 


( \ +i-^v 

27r2 7 o s Wsinh (s) ds^y 

(60) 


The asymptotic behavior for small coupling constant g of the integral is reproduced exactly at the lowest order the 

2 o 

one-loop Effective action. Furthermore, the (unique) non trivial minimum for (1601) is found to be at ^^ 4 ° fa 0.361, as 
it can be read out from Fig.l. 



g‘^9^/k^ 


FIG. 1: Plot of the function (1601) (blue-continuous) and the one-loop (1141) (red-dashed). Notice the presence of a non zero global 
minimum for ^Sr ~ 0.361. 


C. FRG improved YMC Lagrangian 


In order to check whether the effective Lagrangian calculated in (l60l) and derived for a SU(2) YM theory can actually 
explain dark energy, we need to review if our example satisfies the properties we discussed in Section Pi Al 


1 . 


From Fig.l it’s evident that the function (1501) has a non zero global minimum. The exact position of this 

2 0 

minimum can be computed numerically, and in terms of dimensionless quantities is found to be ^^, 4 ° ~ 0.361. 


2. It is possible to compute the asymptotic expansion of Wfe(0) for small values of the YM coupling constant g. 
In this limit we are able to reproduce the one-loop result derived by Savvidy [ 2 ^: if we perform the Taylor 
expansion of 1 / sinh^(s) = | -b O (s^), we find indeed 


yVk{o) = 




1 11 

2 247r2 


g'^OLog 



(61) 


3. yVk{0) shows an asymptotic behavior at least linear in 6. This means that for 9 ^ 
integral tends to 1 and the only 9 dependence is the overall one, namely 


the exponential in the 


Log {yVk{9)) = Log ( 6 >) -b 0{9). 

The condition —?► 1 then follows immediately. 

In the following sections, we show in detail that the YMC described by our toy model evolves from a radiation like 
component to a dark energy one. 
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1. A non-interacting YMC model 

We have already shown that in the case of a non-interacting YMC model the condensate evolution equation is implicitly 
solvable, and that the solution 0{a) can be obtained by inverting the equation 

^f6Wk (0) = aa-^ , (62) 

where a is a coefficient of proportionality that depends on the initial conditions. We can then fix the renormalization 
scale k by comparing the “predicted” YMC fractional energy density with the measured Dark Energy fractional 
energy density (Da = 0.735) finding for a big range of initial conditions k ^ 3.2/i^/^10“^eE. This energy scale, as 
already noted in is low compared to typical energy scales in particle physics, such as the QCD and the weak- 

electromagnetic unification, and this prevents the identification of this YMC as a condensate of some SM gauge fields. 
Then we can study the evolution of the YMC energy density and its EOS for different values of a = (10“^, 10“^), 
still obtaining the same asymptotic values. Results are summarized in Fig. 2. 




Log (2 + 1 ) Log 0-1-1) 

FIG. 2: In the free YMC dark energy model, the evolution of the YMC fractional energy density (left panels) and EOS (right 
panels) for the models with different initial conditions (a = 10“^ on the left and a = 10“® on the right). 


2. Q oc pymc 

For YMC models enjoying an interaction proportional to the YMC energy density, we have already discussed in 
Section HID II the equation defining the existence of a fixed point — derived from the differential equation evolution 
system (P)) - ([57|) — and the condition to be imposed on the coupling parameter a in order to characterize an attractor 
solution. Here we report numerical computations on the position of the fixed point and on the value of the critical 
exponents at the fixed point. The fixed point exists for every positive value of the coupling parameter and is always 
attractive, as shown in Fig. 3. 


IV. CONCLUDING REMARKS 

The query whether the YMC may actually provide a consistent and physically reliable model for dark energy, along 
the lines of the analysis first developed in Refs [IMl, has been addressed in this paper within the attempt of finding 
support to this theoretical hypothesis in the non-perturbative approach to the calculation of the Yang-Mills effective 
action. 

In particular, we have first discussed the properties that the effective Lagrangian yV{6) must posses in order to drive 
the Universe towards a cosmological dark energy phase. If a condensate exists, i.e. if the non-perturbative effective 
action has a minimum in the YMC order-parameter 0, then the model can actually reproduce the dark energy behavior 
of the expanding universe at small redshift. If the effective action scales at least like the bare Yang-Mills action for 
high-energy scale, at high redshifts it entails the EOS of radiation. Moreover, internal consistency also requires that 
perturbative one-loop results must be still recovered in the appropriate asymptotic limit. 

We have then focused on the particular example provided of the SU(2) Yang-Mills bare action. We have deployed 
non-perturbative techniques mutated from the FRG method, in order to show that for the SU(2) Yang-Mills a non¬ 
trivial minimum indeed exists, and that the high-energy scale regime approaches known perturbative results, and 
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^C) 



FIG. 3: For the coupled YMC dark energy models, the plot against the coupling parameter a of the critical order parameter 
and matter density 6c, Xc (top panels) and of the critical exponents Ai,A 2 at the fixed point (bottom panels). In the bottom 
panel, on the left is reported the case in which Q oc Pymc, and on the right the case in which Q oc pvMC + Pm. 


yields radiation dominated EOS. Despite our conclusions can be only based so far on this particular example, this 
successful check of the requirements necessary to have a viable YMC dark energy model seems to us extremely 
encouraging in strengthening the possibility that YMC can work as a model for dark energy. 

The improvement of the non-perturbative techniques may allow in the future to extend the present analysis to 
the cases of S\J{N) gauge-groups, or more in general to other classes of Lie-groups. At this purpose, we may either 
ask ourselves whether condensation can work for any SU(A^) group, or whether consistency of the model necessarily 
predicts a maximal values of N in order the mechanism to work. Not disconnected to these questions, it comes the 
query on the relation between the Yang-Mills fields involved, which are necessary in order the condensate to form, 
and the elementary-particle fields advocated to explain dark matter. Indeed, it would be tempting to try to link YMC 
dark energy models to other models for dark matter, such as the ones referred to in the literature |3ll4^ as Mirror 
Standard Model for dark matter. 
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Appendix A: Action and variations 

We consider the bare action of the Yang-Mills theory 

s = ^l dxFrp;,. 

We split the field into a background plus a fluctuation = A^ -|- a^. The quadratic part of this action in the 
fluctuation field is 
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dx [o^alD^al - D Ucala^) 

= ^jdxal ((5,f,n + D^D’^S.b + gF^'^’^fabc) , 

where the bar quantities are made out of the background fields. To compute the inverse propagator we need to add 
also a gauge fixing action and the corresponding Faddeev-Popov ghost action: 

Sgf = ^ y dxZ)^a“Dya“ , Sgh = j dxDgC^D^c’' 

The FRGE splits into the trace over the transverse part and the longitudinal part of the connection field and the 
ghost sector: 

isTr (^(2) ^ = ^Trr + Rk) dtRk + ^Tr^ (^(2) ^ ^ g^j^^ 

- It , Irp adtRk _ dtRk 

“ 2 + i?fc 2 + aRk □ + i?fc ' 

Calculations are simplified by the fact that in the Landau gauge a —>■ 0 the longitudinal trace can be dropped out 
completely. 


Appendix B: Self-dual field configuration 

We may choose a background field configuration that allows to project onto the effective potential W{9). Hence a 
covariantly constant field strength with = 0 suffices. Since the spectrum of the Laplace-type operators, like 

Dl^'^ = DScbS^'^ + gF°^^'' fabc, or at least the heat-kernel trace for these operators has to be known, we have a limited 
choice in the possible background field configurations. To avoid problems with tachyonic modes, which indicate the 
instability of a background, we project onto the only known stable covariantly constant background, which is self-dual, 
namely F^y = = F^^. Then we set F^i, = 0. Apart from Foi = F 23 = B = const, all other non-zero 

components follow from the antisymmetry of the field strength tensor. Due to the enhanced symmetry properties 
connected to the self-duality, zero modes exist. These carry important information and have to be regularized carefully, 
since the standard choice for Rk is zero on the zero mode subspace. 


spec = 2gBi ( n -f TO -|- 2) with n, TO G N and with multiplicity 2 in 4 dimensions, 

= 2gBi (n -I- to) with n, to £ N and with multiplicity 2 in 4 dimensions, 
spec (□) = 2gBi (n -I- to -|- 1) with n, to G N 


with a degeneracy factor . Herein Bi = \vi\B and vi is given by vi = spec |(T“ I = l| with the 

generators of the adjoint representation T“ and therefore for a general gauge group depends on the direction of the 
unit vector n. 
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